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O ■ ABSTRACT 

Spectral fits to X-ray data from both NGC 253 and M 82 provide ambiguous results. The 
Qh' so-called "best fit" results depend on the instrument with which the data were obtained and 

Q . obviously on the choice of spectral model composition. We show that different spectral models 

\->i \ can be fit equally well to BeppoSAX data of both galaxies. Metallicities are unreliable in general, 

. with a strong dependence on the choice of model. Preference to one particular spectral model 

can only be given by combining spectroscopic and imaging X-ray data from all available satellites 
{ROSAT, ASCA, and BeppoSAX). Based on spectra of NGC 253, we demonstrate that a model 
consisting of two or more thermal plasma components plus a hard power law continuum and Fe 
Ka line emission can explain all observations. These model components represent the integral 
' spectrum of thermal gas and compact sources in starburst galaxies that are most likely supernova 

remnants and X-ray binaries. The same model can fit the X-ray data of M82, but there the 
evidence, from i?05'j4T imaging, of the existence of compact sources which might represent high- 
mass X-ray binaries is weaker. This implies that its hard X-ray emission, which is extended in 
i? OS'^T images, might-if truly diffuse-be dominated by a very hot (several keV energy) thermal 
gas component. 

Subject headings: Galaxies: individual: (NGC 253, NGC 3034 — M82) — Galaxies: halos — Galaxies: 
starburst — Galaxies: ISM — X-rays: galaxies 

1. Introduction land et al. 1997; Tsuru et al. 1997; Vogler & 

Pietsch 1999). All authors agree on the general 
There have been a number of papers recently complexity of the X-ray properties, but-depending 
discussing ASCA and ROSAT X-vay observations the data and on the spectral models used for 
of nearby starburst galaxies, in particular the two ^.^e spectral analysis-they reached different con- 
best-studied systems NGC 253 and M 82, with dif- elusions. Especially the choice of spectral model 
ferent interpretations of the resulting spectral fits components and the resulting best-fitting element 
(Moran & Lehnert 1997; Ptak et al. 1997; Strick- 
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abundances are under debate. 

By treating all available imaging and spec- 
troscopy data from ASCA and ROSAT (both HRI 
and PSPC) in a self-consistent manner, Dahlem 
et al. (1998; hereafter Paper I) and Weaver et al. 
(2000; hereafter Paper II) were able to reconcile 
the apparent discrepancies based on a mini-survey 
of 5 nearby edgc-on starburst galaxies. 

The results from Paper I and Paper II suggest 
that the combined ASCA and ROSATPSPC inte- 
gral spectra of NGC 253 and M 82 can be fit with 
comparable values of by different combinations 
of spectral components, which means that there is 
an ambiguity in the choice of the best-fitting spec- 
tral model. By cross-checking the spectral results 
with ROSATPSPC and HRI imaging data, a spec- 
tral composition of (at least) two thermal plasmas, 
with temperatures in the ranges 0.1-0.4 keV and 
0.6-0.8 keV, respectively, plus a hard power law 
component turns out to be the only model combi- 
nation that can explain all observational data of all 
galaxies in the sample simultaneously (Paper II). 

This is in contrast with the recent findings by 
Cappi et al. (1999; hereafter C99), based on Bep- 
poSAX data of NGC 253 and M 82. These authors 
claim that there is "compelling evidence for the 
presence of an extended hot thermal gas" of sev- 
eral keV temperature in these two galaxies. The 
purpose of the current letter is to investigate this 
apparent discrepancy between their results and 
ours by re-analyzing the BeppoSAX observations 
of NGC 253 and M82, taking into account earlier 
results based on ROSAT and ASCA observations. 

2. Observations and data reduction 

All parameters of the BeppoSAX observations 
of M82 and NGC 253 are as described by C99. 
The data were reduced in the standard fashion, 
using SAXDAS 2.0. LEGS and MEGS data were 
extracted for joint spectral fitting from a circu- 
lar region centered on the position of the sources, 
using radii of 4'. Background subtraction was per- 
formed using standard files (Parmar, Oosterbroek, 
& Orr 1999 ***check***). 

Spectral fitting was performed using XSPEC in 
the following way. We first used the input model 
preferred by G99 to ensure that we can reproduce 
their results. Then we tried the model used by us 
for the analysis of the joint ASCA+ROSATPSPC 



spectrum (Paper I and Paper II). 

In the following we will list and discuss our re- 
sults for both NGC 253 and M82. However, since 
we used the same data extraction, reduction, and 
spectral fitting technique for both galaxies, only 
one (NGC 253) will be presented in figures. 

3. Results and discussion 
3.1. NGC 253 

We could reproduce the results by C99 within 
the uncertainties, using a spectral model with 2 
Mekal plasma components (2M; see their Fig. 4). 
The fit to the BeppoSAX data of NGC 253 follow- 
ing the model preferred by us, with a Mekal and 
a power law (M-l-P) component, is displayed in 
Fig. 1. This is evidently also an acceptable fit. 
The goodness of fit for our preferred model (Pa- 
per II) is = 261.9 for 265 d.o.f. {x^/v = 0.99), 
while that for the 2 Mekal (hereafter "2M") model 
favored by C99 is = 282.5 for 268 degrees of 
freedom (d.o.f.) and thus x^/v = 1-05. The M-l-P 
model fits the data better than the 2M model at 
the highest and lowest energies of the passband. 
The results of the two spectral fits to the Bep- 
poSAX data of NGC 253 are tabulated in Table 1. 
All uncertainties are given at the 90% confidence 
level for one interesting parameter; note that these 
only apply under the assumption that the choice 
of model components represents the different con- 
tributing emission mechanisms correctly. 

The softest thermal emission component found 
in the -ROS'^Tdata is not required. When includ- 
ing a thermal plasma of 0.26 keV energy, x^ is 
improved, but not significantly. Therefore it was 
left out in the fits to the BeppoSAX data. 

The hard part of the X-ray spectrum can be fit 
with a power law that is compatible with those 
of Galactic X-ray binaries (XRBs) and can thus 
be explained naturally as the continuum emission 
from HMXRBs (Paper II). The integral spectrum 
of all the point sources from the ROSAT PSPC 
observations is consistent with this interpretation, 
as is the contribution of this spectral component 
to the total X-ray flux (Paper I). There is no rea- 
son why all compact sources should emit a thermal 
spectrum. 

The claim that there is a hot thermal gas at a 
few keV energy (C99) hinges only on the assump- 
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tion that this is the only mechanism that could 
explain the observed Fe line around 6.7 keV. How- 
ever, supernova remnants (SNRs) and XRBs, in- 
cluding high-mass XRBs (HMXRBs), can produce 
both fluorescent and thermal Fe line emission, i.e., 
at 6.4 keV and 6.7 keV, respectively (e.g., Nagase 
1989; White, Nagase, & Parmar 1995; Licdahl et 
al. 1999), while C99 assume that all line emission 
is of thermal origin. Thus, the observed Fe line 
might well be a superposition of emission from hot 
gas and from X-ray binaries. The fitted equiv- 
alent width of the Fe line (Table 1) is in agree- 
ment with this interpretation. It, too, might be a 
composition of a (broad) thermal component and 
a (narrow) XRB contribution. However, a com- 
posite line fit cannot be performed based on the 
current data, because the line is only marginally 
resolved. 

The above finding that different models can fit 
the data equally well proves that the "optimal" fit 
is model-dependent, as already stated in Paper II. 
Therefore, there is no reason to reject the M-l-P 
model. Moreover, it is, as argued in Paper II, the 
physically most plausible model choice. 

Taking into account the ROSAT imaging re- 
sults, which indicate clearly that there is a consid- 
erable number of unresolved compact sources in 
the central part of the disk of NGC 253 (Paper I), 
the most likely identification of these sources- 
based on their spectral properties and soft X-ray 
luminosities, ix-is that they constitute a popula- 
tion of HMXRBs (Paper I). Thus, part of the emis- 
sion distribution seen by BeppoSAX is not truly 
"diffuse", but smeared out by its broad point- 
spread function. 

These point sources detected by ROSAT con- 
tribute about 50% of the flux from the central 
disk (Paper I and Paper II). Thus, they are very 
significant contributors to the measured total flux, 
especially in the hard part of the X-ray spectrum. 

On the other hand, the spectral model preferred 
by C99 does not take into account the presence 
of HMXRBs and their spectral signature. Given 
the luminosities of emission mechanisms tracing 
the presence of high-mass stars in galaxies like 
NGC 253 and M82, especially far-infrared radi- 
ation, a large number of HMXRBs must be ex- 
pected to be present in them. 

It is still unclear how the previously detected X- 



ray emitting thermal plasma (with temperatures 
in the range of a few tenths of a keV) is heated, es- 
pecially in the galaxy halos, up to several kpc away 
from the disk planes of the starbursts. The pres- 
ence of another, extremely hot medium of several 
keV energy contributing of order 2/3 of the total 
2 10 keV flux, as suggested by C99, would further 
excruciate the problem of energy supply. 

When taking into account the trade-off between 
metallicities and absorbing HI column densities 
in fitting the softest part of X-ray spectra (which 
cannot be resolved by BeppoSAX data only, but 
requires the low-energy response of ROSAT), ex- 
treme subsolar metallicities, Z, are not required 
to obtain a good fit (Paper II) . This iVn vs. Z di- 
chotomy is another, independent ambiguity in the 
minimum space of the spectral fits. Low metal- 
licities in starburst galaxies, i.e., the galaxies with 
the highest star formation rates in the local Uni- 
verse, would be hard to understand because of the 
proven presence of large numbers of massive stars, 
which are the most prolific producers of metals. 

3.2. M82 

The same ambiguities are present in fits to the 
BeppoSAX data of M82. We could fit almost 
equally well the model by C99 and ours from Pa- 
per I and Paper II. Just as for NGC 253, the M-l-P 
model fits the data points at the very highest en- 
ergies slightly better than a 2M model. The good- 
ness of fit is 517.4/446 d.o.f. = 1.16 (2M model) 
and 466.9/442 d.o.f. = 1.06 (M+P model), respec- 
tively. Note that, just as for the combined ROSAT 
+ ASCA data, the BeppoSAX data require an- 
other, soft thermal component to be added to the 
M-l-P model. 

With the M-l-P spectral model composition, we 
obtain almost equally good fits with two very dif- 
ferent metallicities. In one case, Z = 17 Zq (con- 
strained to be > 2.2 at the 90% confidence level), 
in the other Z = 0.13^0. In the high-metallicity 
case the flux at ~ 1 keV is modeled primarily as Ne 
and Fe-L line emission, while in the low-metallicity 
case it is modeled as a peak in the thermal dis- 
tribution. The energy resolution of the LEGS of 
~ 200 eV (FWHM) at 1 keV is insufficient to 
discriminate betwecin the two options. Note that 
these two fits do not yet take into account the 
additional information obtained with ASCA and 
i? 05^7 requiring an additional soft thermal com- 
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ponent (Paper I and Paper II). 

There is less evidence from i? 05^ T imaging for 
the existence of large numbers of compact sources 
in M82. Instead, there appears to be a spa- 
tially extended, hard spectral component. Part 
of this might be truly diffuse, in which case the 
most likely interpretation is that of a very hot 
gaseous component, as suggested by C99. Only re- 
cently the Chandra image by Griffiths et al. (2000) 
showed that there is indeed a population of com- 
pact sources in M82, surrounded by diffuse emis- 
sion. The compact sources in the central part of 
M82 could not be resolved by ROSAT, because 
they are too close to each other. The most likely 
identification is again that they are HMXRBs 
(Griffiths et al. 2000). Individual HMXRBs could 
also explain the observed X-ray variability in the 
hard part of the spectrum, while there is no evi- 
dence in the Chandra data for the presence of an 
AGN (Ptak & Griffiths 1999, Matsumoto & Tsuru 
1999, Gruber & Rephaefi 1999, C99). 

The measured position and equivalent width of 
the Fe fine in M82 of 6.63 ± 0.21 keV and 60 ± 
40 eV, respectively, leaves open whether the line 
emission comes from either binaries or diffuse hot 
gas or a superposition of both. In both M 82 the 
width of the Fe line near 6.6 keV is unresolved. 
Thus, except for the (poorly constrained) position 
of the line centroid, no further information on the 
relative contribution of thermal or fluorescent line 
emission can be made based on the existing data. 
Both model compositions tested above fit the data 
(statistically) so well that no useful constraint can 
currently be derived on the possible contribution 
of both a hot thermal plasma and HMXRBs to the 
2 10 kcV flux of M 82, when fltted simultaneously. 

4. Summary 

There are several ambiguities in the flts to com- 
plex X-ray spectra of starburst galaxies, such as 
NGC253 and M82. The "best-fitting" model is 
not necessarily unique, because the spectral mod- 
els required to explain all observations are more 
complex than can be fit unambiguously to one sin- 
gle dataset. In such cases statements that a fit is 
good at a certain significance level can be mislead- 
ing, because they only apply if the correct spectral 
model composition was chosen. There are also in- 
trinsic degeneracies, i.e., trade-offs of different fit 



parameters against each other (e.g., TVh vs. Z). 

This study makes it clear how important it 
is to consider all available information, including 
in particular X-ray imaging results of extended 
sources, when interpreting their integral spectral 
properties. The new generation of X-ray satellites, 
Chandra, XMM, and Astro- E, will resolve much of 
this ambiguity because of their high spectral res- 
olution, combined with high sensitivity and good 
imaging capabilities over wide bandpasses, render- 
ing possible spatially resolved spectroscopy of indi- 
vidual (classes of) sources within nearby galaxies. 
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Table 1 

Comparison of fit results to BeppoSAX data of NGC 253 



Source^ 


ModeP 


(keV) 


Nn.m 
(10^2 




kTh/r 
(keV/...) 


Nn.h 
(1022 cm-2) 


Zh 

(Ze) 




C99 

Paper II 


2M 
M+P 


0.75 ±0.06 


fi OO+0.22 
0.013(f) 


0.25t°;°^ 
0.16 ±0.02 


5.90t°-^ 
2.22 ± 0.04" 


n OO+0.22 

1.11±0.11 


0.25tH^ 


1.05 
0.99 



^C99 = Cappi et al. (1999); Paper II = Weaver et al. (2000) 
22M = Mekal + Mekal; M+P = Mekal + Power Law 

■^Subscript m for "medium" energy thermal component (Paper II), subscript h denotes the alleged "hot" thermal 

component (C99) 

"Fc Ka line fitted independently at 6.62 ± 0.07 keV, with an equivalent width of EW = 480 ± 110 eV 
^Parameter fixed 
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